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Abstract: The analysis of data from the MODIS instrument 
on the Terra platform to derive global distribution of aerosols 
assumes a set of relationships between the blue, p blue , the red, 
p red , and 2. 1 micrometers, p 21 , spectral channels. These 
relations have been established from a series of measurements 
indicating that p blue ~ 0.5 p red ~ 0.25 p 2J . Here we use a model 
to describe the transfer of radiation through a vegetation canopy 
composed of randomly oriented leaves to assess the theoretical 
foundations for these relationships. The influence of varying 
fractional vegetation coverage is simulated simply as a linear 
combination of pure soil and pure vegetation conditions, also 
known as Independent Pixel Approximation (IP A). Calculations 
for a wide range of leaf area indices and vegetation fractions 
show that p blue is consistently about 1/4 of p 2 A as used by 
MODIS for the whole range of analyzed cases, except for very 
dark soils, such as those found in bum scars. For its part, the 
ratio p r ed/p 2 .i varies from less than the empirically derived value 
of 1/2 for dense and dark vegetation (p 2 j<0.1), to more than 1/2 
for bright mixture of soil and vegetation. This is in agreement 
with measurements over uniform dense vegetation, but not with 
measurements over mixed dark scenes. In the later case, the 
discrepancy is probably mitigated by shadows due to uneven 
canopy and terrain on a large scale. It is concluded that the value 
of this ratio should ideally be made dependent on the land cover 
type in the operational processing of MODIS data, especially 
over dense forests. 

1. Introduction 

A new generation of space-borne Earth observation 
instruments has been launched and can be used to monitor 
aerosols over land. Satellite sensors differ in their measurement 
capability: MODIS on the Terra satellite is using a wide spectral 
range (0.47 to 2.1 pm) to derive the aerosol optical thickness 
and properties [Kaufman et al , 1997a]; MISR, also on Terra, 
retrieves aerosol properties using a modest spectral range (0.44 
to 0.86 pm) but observing the Earth from nine different 
directions simultaneously along the satellite track [Diner et al, 
1989, 1998]. The POLDER instrument, operating in 1996-1997, 
using the same spectral range, observed the Earth from a wide 
range of angles in two dimensions and derived the aerosol 
properties with the help of polarization. The ATSR instrument 
combines 4 spectral channels (0.55 to 1.65 pm) with two view 
directions [Veefkind, 1998]. In all these instruments the 
measured signal consists of sunlight reflected by both the Earth 
surface and the aerosol layer. Some assumptions about the 
surface reflectance properties and aerosol properties have 
therefore to be made to separate the aerosol contribution to the 
signal from that of the surface. Each instrument takes advantage 
of its unique configuration to derive useful information on 
atmospheric aerosols. 


In this paper, we concentrate on the MODIS technique 
based on the assumption that the ratios of the surface 
reflectances in the blue (0.47 pm) and red (0.67 pm) to the 
reflectance at 2. 1 pm are constant in space and time at the fixed 
values of 0.25 and 0.50, respectively. Our simulations were 
carried out at 0.44 pm, which has similar characteristics to the 
MODIS 0.47 pm channel. This ratio was established using 
remote sensing measurements in the United States [Kaufman et 
al. , 1997b] and confirmed in an independent experiment in Israel 
[Karnieli et al., 2001]. The global application of the technique 
[Chu et al., 2001] shows that the method works well in most 
cases. However no theoretical basis for this assumption had been 
developed so far. More specifically, the MODIS technique is 
based on the fact that sub-micrometer aerosol particles are 
transparent at 2.1 pm (i.e., smoke and pollution but not dust). 
The MODIS algorithm first considers the reflectance at 2. 1 pm 
at the top of the atmosphere, and then uses this reflectance value 
to estimate the surface reflectance in the blue and red channels, 
using the empirical relationships described above. The aerosol 
optical thickness and properties are derived directly from the 
measured excess reflectance in the blue and red channels. In this 
paper, we provide some theoretical basis in support of the fixed 
spectral ratio approach, using surface bi-directional reflectance 
model of Gobron et al. [1997]. 

2 Model simulations 

We simulate the angular and spectral properties of fully and 
partially vegetated surface using the plant canopy radiation 
transfer model developed by Gobron et al, [1997]. This model 
solves the radiation transfer problem in the case of finite-size 
oriented leaves uniformly distributed over a scattering soil. The 
leaves are not clumped into trees, but are distributed in the space 
as illustrated in Fig. 1. The model computes the Bidirectional 
Reflectance Factors (BRFs) for mixed surfaces using a linear 
mixing approach: 

Pmixed(^o>^) = ^|fiPi(^o>^) (1) 

Here, fj is the relative fraction of the pixel area occupied by 
surface type i. Allowed surface types include either vegetation 
with a given leaf area index (LAI) or soils, while £2 0 and 
denote the solar illumination and satellite observation directions, 
respectively. Table 1 summarizes the leaf and soil spectral 
properties. The latter were extracted from the database of Price 
[1995] and represent a wide range of actual soils. The leaf 
spectral properties are based on the leaf model of Jacquemoud 
andBaret [1990] for typical standard conditions of chlorophyll 
concentration and water content in a healthy leaf. The radiative 
transfer model has been extensively evaluated and tested against 
a set of three-dimensional models [Pinty et al., 2001]. It was 
shown to be one of the best performing models within the 
structurally homogeneous model category. The linear 



assumption used in eq. 1 to represent the angular reflectance of 
mixed pixels has been proven to be quite accurate for medium to 
low spatial resolution sensors. The calculations are performed 
for a wide range of leaf area index (0.05 to 5), vegetation 
fraction (0.25 to 1), view zenith angle (0° to 40°) and solar 
zenith angle (1° to 51°). The relative azimuth between the Sun 
and the observer is 0°, 90° or 180°. 



Figure 1: Artist view of the structurally homogeneous leaf 
distribution implemented by the canopy radiation transfer model 
(top), and of the mixed surfaces used in our simulations, here for 
a fractional vegetation cover of 50% (bottom). 


Table 1: Spectral properties of leaves and soils used in the 
simulations 


Reflective surface / wavelength 

0.44 pm 

0.67 pm 

2.1 pm 

Leaf reflectance 

0.054 

0.054 

0.148 

Leaf transmittance 

0.005 

0.014 

0.172 

Dark soil reflectance 

0.068 

0.126 

0.140 

Medium soil reflectance 

0.102 

0.238 

0.357 

Bright soil reflectance 

0.139 

0.351 

0.536 


3. Results 

The plant canopy radiation transfer model [Gobron et al , 
1997] is used to calculate the spectral reflectance of mixed 
scenes using eq. 1, where f v represents the fractional vegetation 
cover (f v = 0.25, 0.5, 0.75 and 1.0) and (1- f v ) is the fractional 
soil cover. The leaf area index (LAI) varies from 0.05 to 5. 
Throughout the ensemble of simulations, the vegetation density 
(f v *LAI) thus varies from 0.01 to 5. Note that shadows due to 
variability of the height of vegetation and terrain relief are not 
accounted for in these calculations. We shall evaluate the 
implications of this issue in the discussion section. Figure 2 
shows scatter plots of the simulated mixed surface reflectance at 
0.44 (iun and 0.67 pm as a function of the reflectance at 2. 1 pm 
for all cases corresponding to medium and bright soils. The 
vegetation density is represented in color, ranging from dark 
green for f v *LAI=5 to red for f v *LAI=0.01. Lines corresponding 
to the ratios of 1/4 and 1/2, for the blue and red respectively, as 
assumed in the MODIS algorithm [Kaufman et al, 1997a] are 


superimposed. Quasi-linear relationships are observed between 
the blue (0.44 pm) or red (0.67 pm) and the 2.1 pm channels for 
a wide range of simulated conditions. The simulated 0.44 to 2. 1 
pm ratio closely follows the assumed ratio of 1/4 over the whole 
range of surface reflectance values. In the case of the 0.67 over 
2. 1 pm ratio, however, the simulated ratios are lower than 1/2 
for dark surfaces and larger than 1/2 for bright surfaces. The 
figure shows a transition from ratios p re d/p 2 .i around 0.25 for 
dense dark vegetation (typically forest), to ratios larger than 1/2 
for bright soils (see Table 1). In fact, measurements do show 
that for dense homogeneous forests the relationship between the 
reflectance in the 0.44, 0.67 and 2. 1 pm is such that p blue ~ p red ~ 
0.5p 2 .i [Remer et al, 2001; Gatebe et al, 2001], as predicted by 
the model. For mixed dark scenes, however the measurements 
show that p blue ~ 0.5 p red ~ 0.25 p 2 .i [Kaufman et al, 1997a; 
Karnieli etal, 2001]. 



Figure 2: Scatter plots of the simulated mixed surface 
reflectance at 0.44 pm (top) and 0.67 pm (bottom) as a function 
of the reflectance at 2.1 pm for medium and bright soils. The 
vegetation density is represented by the color bar, varying from 
dark green for f v *LAI=5 to red for f v *LAI=0.01. f v is the 
vegetation fraction, and LAI is the leaf area index. Lines for the 
ratio of 1/4 and 1/2, assumed in the MODIS algorithm [Kaufman 
etal, 1997a] are also shown. 


Figure 3 shows similar calculations to those of Fig. 2, but 
this time including simulations with dark soils. The plots are for 
a limited reflectance range of 0 to 0.2 to emphasize the range 
that is used by the MODIS analysis. Note also that the color bar 
indicates the various LAI values instead of the vegetation 
density (f v *LAI). As already shown in Table 1, dark soils reduce 
the reflectance at 2.1 pm much more than in the red or blue 
wavelengths, thus resulting in a much steeper slope. Such dark 
soil conditions are typical, for instance, of recent burn scars. 
Other soil conditions intermediate between dark and medium, 
would probably end up with a series of different slopes. 



Surface reflectance at 2*10 


Figure 3: Similar calculations to those in Fig. 2, but including 
the dark soil reflectance (dots) from Table 1. The scale of the 
surface reflectance is limited here to a lower range of 0 to 0.2. 
The colors represent the leaf area index, LAI, which varies from 
0.05 to 5. For each LAI value, the vegetation fraction varies 
from 0.25 to 1.0. 

It is interesting to see whether the ratios of Fig 2 depend on 
the view and solar illumination directions. The ratios of the blue 
and red reflectances to the 2. 1 pm reflectance are shown in Fig. 
4 for two solar zenith angles: 21° and 51°, representing 
conditions of high and low Sun. They are represented as a 
function of the view zenith angle. The figures correspond to full 
and half fractional vegetation covers, for LAI varying from 0.05 
to 5 and for conditions where the reflectance at 2. 1 pm is less 
than 0.25. It can be seen that the ratios are weakly dependent on 


the view zenith angle for low solar zenith angle. They are also 
weakly dependent on view zenith angle for vegetation mixed 
with soil. The dependence is stronger for higher solar zenith 
angles over full vegetation cover. MODIS will encounter high 
solar zenith angles at high latitudes, e.g., over the boreal forests. 



50% Coverage- 
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Figure 4: Ratios of the blue (0.44 pm) and red (0.67 pm) 
reflectances to the 2. 1 pm reflectance, represented as a function 
of the view zenith angle for solar zenith angles of 21° and 51°, 
typical of high and low Sun. The figures are given for full (top) 
and half fractional vegetation cover (bottom), for conditions 
where the reflectance at 2. 1 pm is less than 0.25. The LAI varies 
from 0.05 to 5 following the color convention used in Fig 3. 

4. Discussion and conclusions 

A model describing the transfer of radiation through a 
vegetation canopy composed of randomly oriented leaves is 
used to formulate a theoretical basis for the empirical 
relationship between the surface reflectances in the blue, red and 






2. 1 pm channels assumed in the MODIS aerosol algorithm. The 
calculations were carried out for a wide range of leaf area index 
and vegetation fractions, using a linear mixing between 
vegetation and soil reflectances. 


We show that the reflectance in the blue is consistently 1/4 
of the reflectance at 2. 1 pm for all analyzed cases over medium 
and bright soils. This finding extends the limited empirical basis 
of the MODIS aerosol algorithm to many environments found 
around the globe. 

The ratio of the reflectance in the red to that at 2.1 pm 
varies from slightly less than the empirical value of 1/2 for 
dense and dark vegetation (p 21 <0.1), to slightly more than 1/2 
for situations involving vegetation mixed with bright soils. This 
is in agreement with aircraft measurements over uniform dense 
vegetation [Remer et al, 2001, Gatebe et al, 2001] and over 
exposed soils but not for cases where vegetation is mixed with 
dark soils. One possible explanation for the differences for 
mixed scenes is the effect of shadows in an uneven terrain or 
uneven height of vegetation. Neither effect is explicitly included 
in the present model calculations. For example, a reflectance at 
2.1 pm of p2.i=0.1 corresponds to a surface reflectance of 
p red =0.035 in Fig. 2 or 3. However using values from Fig. 2 we 
can construct a surface cover with a combination of 60% of dark 
surface (p 2 .i=0.06 and p red =0.015) with 40% of bright surface 
(P2. i=0.30 and p red =0.18). If this combination is 35% shadowed 
and if these shadows are totally black in both channels, the 
linearly composited reflectance gives the same reflectance at 2. 1 
(P2.i=0-10) but higher reflectance in the red of p red =0.052, than 
the original value of p red =0.035. 

Nevertheless, and in agreement with the mentioned aircraft 
measurements, the conclusion of this study is that for very dark 
surfaces (e.g., p 2 .i < 0.06) the surface reflectance in the red 
channel, may be overestimated by assuming a ratio of 1/2, 
resulting in an underestimate in the derived aerosol optical 
thickness. A smaller ratio should be used however more 
extensive field and modeling experiments have to be conducted 
in order to better quantify this ratio. 
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Popular Summary 

We derive routinely the concentration of aerosol particles, their properties and 
effect on solar radiation from satellite data, namely the MODIS instrument on the 
Terra platform. Aerosol is the haze part of urban pollution, smoke forest fires or 
desert dust. The aerosol information is derived over ocean and land. While ocean 
is uniformly dark and thus easy background to see the reflection of sunlight from 
the aerosol, the problem is more complex over the land. MODIS detects sunlight 
reflected both by the land surface and by the aerosol. How to distinguish 
between them? 

We found that pollution and smoke aerosol is transparent at long invisible solar 
wavelengths. We therefore can use them to see how what the surface is 
composed of without aerosol interference, and try to use this information to 
guess how reflective the surface is in the visible part of the spectrum. Inferring 
the surface reflection of sunlight in the visible from the long wavelengths is 
based on collection of many measurements in different parts of the world that we 
did before the launch of MODIS. However we never had a clear physical 
understanding of this method. In this paper we for the first time use a model that 
simulates surface reflection across the solar spectrum for mixtures of soils and 
vegetation. We prove that our relationship between the surface properties in the 
visible and in the long wavelengths works most of the time, and found methods 
to correct them in cases it does not work. 



